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ABSTRACT 

Due  to  increasing  application  of  Micro-Electro-Mechanical  Systems  (MEMS),  are  in  constant  need  of  better 
and  improved  understanding  of  fluid  flow  and  heat  transfer  characteristic sf  in  micro  channels.  This  paper  offers  a 
computational  analysis ,  by  using  the  FLUENT  14.0  and  ICEM  software.  The  purpose  of  this  study  is  to  explore ,  the 
Thermo -hydraulic  traits  of  single -phase  fluid  flow  by  using  air ,  liquid  nitrogen ,  and  liquid  hydrogen  as  flowing  fluids  in 
smooth  circular  micro  channels  over  the  hydraulic  diameter  range  from  1  pm  to  100  pm.  The  mass  flow  rate  varies 
from  0.001  kg/sec  to  0.01  kg/sec  for  every  particular  diameter.  The  length  of  micro  channel  was  kept  constant  as  0.1 
meter.  Flow  in  the  micro  channels  was  considered  as  turbulent  and  the  variation  in  Reynolds  number  value  for  all  three 
flowing  fluids  is  considered  to  primarily  examine  the  diversification  in  friction  factor  and  nusselt  number ,  over  the 
diameter  range.  Leisurely  contemplation  is  made  to  justify  the  parity  of  computational  results  from  conventional  theory , 
while  friction  factor  and  pressure  drop  was  seen  to  be  compatible  with  classical  theory  of fluid  flow  for  all  three  flowing 
fluids. 
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INTRODUCTION 

The  recent  advancement  in  micro  heat  exchanger  technology  has  pushed  the  thermal  restriction  in 
electronic  [4]. 

This  idea  of  implementing  the  microchannels,  for  quenching  the  heat  emitted  from  very  large-scale 
integration  (VLSI)  circuits  was  first  proposed,  by  Tuckerman  [5].  He  expounded  about  the  direct  contact  and  cold 
plate  hybrid  technology,  by  exerting  laminar  flow  of  water  in  the  microchannels  and  imbibing,  150  W/cm2  of  heat 
flux.  The  study  carried  out  [6-8],  to  depict  the  flow  and  heat  transfer  characteristics  of  water  flowing  in  the 
microchannels  of  hydraulic  diameter,  ranging  from  130pm  to  340pm.  It  was  observed  that,  laminar  heat  transfer 
took  place  at  Re=200-700  and  fully  turbulent  convective  heat  transfer  was  achieved,  at  Re=400-1500. 

In  the  experiment,  performed  [9-10]  to  investigate  the  flow  characteristics  in  the  microchannels,  Isopropyl 
data  demonstrated  that,  friction  factor  was  not  dependent  of  Re,  but  it  was  dependent  on  the  size  of  the  channel. 
However,  in  case  of  silicon  oil,  friction  factor  was  independent  of  size  of  the  channels,  but  it  was  merely  dependent 
on  the  value  of  Re.  The  investigation  conducted  [11]  for  the  liquid  flow  in  round  and  square  microchannels  of 
hydraulic  diameters,  ranging  15-150  pm  over  a  Reynolds  number  range  of  8-2300  and  it  was  seen  that,  there  was 
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no  distinguishable  deviation  from  stokes  theory,  which  suggests  that,  the  flow  in  microchannels  strongly  can  be  predicted 
by  the  conventional  theory. 

Study  conducted  on  two  phase  pressure  drop,  in  horizontal  circular  microchannel  [12]  having  hydraulic  diameter 
781pm  showed  that,  Two  phase  pressure  drop  increases  linearly  with  heat  flux  and  vapour  friction,  for  both  R134a  and 
R245fa  refrigerants.  Another  study  [13]  performed  an  Experimental  and  Numerical  Study,  of  Flow  and  Heat  Transfer  in 
Trapezoidal  Microchannels.  It  was  indicated  that,  for  a  particular  pumping  power,  if  heating  power  is  increased  at  the 
substrate,  the  thermal  resistance  decreases,  however  the  sub-thermal  resistance  proportion  is  nearly  unchanged. 

Amirah  et  al.  [23],  investigated  the  effects  of  hydraulic  diameter  and  aspect  ratio,  in  flow  and  heat  transfer  in 
microchannels.  The  range  of  hydraulic  diameters  was  0.1-1  mm  and  the  aspect  ratio  was  kept  1.  Reynolds  number  ranger 
was  100-2000.  Results  of  the  study  suggested  that,  friction  factor  decreases  with  the  aspect  ratio,  greater  than  2  and  also 
found  that;  aspect  ratio  doesn’t  affect  the  heat  transfer  co-efficient.  Another  study  conducted  by  Omer  et  al.  [24],  to 
investigate  the  conjugate  heat  transfer  performance  of  a  heat  sink.  Results  predicted  significant  increase  in  heat  transfer  but 
secondary  flow  was  observed  to  be  blocking  the  heat  transfer.  In  an  another  investigation  conducted  by  Seyed  Ebrahim 
Ghasemi  et  al.  [25],  investigated  the  heat  transfer  performance  of  a  heat  sink,  used  for  the  cooling  of  electronic 
applications.  In  this  experiment,  the  main  channel  of  hydraulic  diameter  4  mm,  6  mm,  8  mm  was  investigated.  Results  of 
the  study  suggested  that,  the  heat  sink  having  Dh=  4  mm  has  much  lower  thermal  resistance,  than  the  rest  and  increase  in 
diameter  decreases  in  pressure  drop. 

TRANSITION  TO  TURBULENCE 

Conventional  theory  suggests  that  when  Reynolds  number  increases  to  a  certain  value  the  fluid  flow  in  duct 
displays  an  astonishing  transition  from  laminar  flow  in  to  turbulent  flow.  The  cause  of  generation  of  turbulence  in  flow  is 
eddies  formation  due  to  high  fluid  velocity.  The  fluid  flow  inside  a  smooth  pipe  having  uniform  cross-section  and  low 
Reynolds  number  results  fluid  particles  to  move  along  the  streamline.  According  to  the  Newton’s  law,  viscous  forces  try  to 
slowdown  the  particle  near  to  the  wall.  The  flow  is  in  well  order  and  called  as  laminar  flow.  However,  at  higher  Reynolds 
number  this  well  ordered  flow  cease  to  exist  and  mixing  of  particle  from  different  layers  starts.  The  fluctuation  in  fluid 
velocity  at  point  causes  momentum  exchange  in  transverse  direction.  The  transition  from  laminar  to  turbulent  is  huge  in 
flow  resistance  pattern. 

The  basic  governing  equations  [14]  are, 

•  Continuity  equation 

PiV-U  =  0  (1) 

•  Momentum  equation 

Pl(.U.AU)=  -VP+  V.Gu,V[/)  (2) 

•  Energy  equation 

plcPil(U.VT)  =  klV2T  (3) 

A  study  performed  on  general  boundary  condition  for  liquid  flow  at  solid  surfaces  [15]  provided  the  molecular 
dynamics  simulations  to  calculate  the  slip-flow  boundary  condition  dependence  on  shear  rate. 


Impact  Factor  (JCC):6.8765 


NAAS  Rating:  3.11 


Computational  Investigation  on  Thermo  hydraulic  Characterization  of  Liquid 
Hydrogen  and  Liquid  Nitrogen  in  Micro  channels 


123 


Aliw  Ufiuid  U-wall  Ls  Qy  (4) 

du 

Where,  Ls  is  fixed  length  of  slip,  and  —  is  the  rate  of  strain  at  the  wall.  The  objective  of  Thompson  and  Troian's 

numerical  analysis  was  to  quantify  the  ascendancy  of  slip  at  a  solid-liquid  interface,  as  the  interfacial  parameters  and  the 
change  in  shear  rate. 

RESULTS  AND  DISCUSSIONS 

Present  2-D  computational  study  on  Thermohydraulic  characterization,  was  performed  to  observe  the  behaviour 
of  fluid  flow  and  heat  transfer  parameters,  in  microchannels.  In  the  investigation,  the  microchannels  of  varying  diameter 
from  1pm  to  100pm,  were  modelled  as  keeping  the  length  0.1  meter  constant.  The  air,  liquid  nitrogen  and  liquid  hydrogen 
were  chosen  as  flowing  fluids.  For  every  diameter,  a  wide  range  of  mass  flow  rate,  from  0.5  gram/sec  to  10  gram/sec  were 
used  for  the  analysis  of  fluidic  and  thermal  characteristics. 

This  study  is  based  on  following  assumptions. 

•  Flow  in  microchannel  is  turbulent. 

•  Slip  at  the  wall  is  zero  (no  slip). 

•  Axi-symmetric  condition. 


Figure  1:  Comparison  of  Dimensional  Entrance  Length  of  Air,  Liquid  Hydrogen  and 
Liquid  Nitrogen  at  Constant  Mass  Flow  Rate 

Hydraulic  Characteristics 
Entrance  Length 

The  graphs  for  the  numerical  dimensional  entrance  length  vs.  Reynolds  number  in  the  figure  1,  was  plotted  for  the 
hydraulic  diameter  range  of  1pm  to  100pm  over  the  constant  mass  flow  rate  of  1  gram/sec  for  the  comparison  of  three 
flowing  fluids  as  air,  liquid  hydrogen  and  liquid  nitrogen  for  circular  geometry  of  the  microchannels.  Fluid  flow  velocity  at 
inlet  was  considered  to  be  uniform.  According  to  the  figure  1,  the  dimensional  entrance  length  does  not  have  significant 
impact  for  the  increasing  value  of  diameter  and  it  was  seen  to  be  getting  independent  of  the  diameter  as  well  for  all  three 
fluids. 

Friction  Factor 

A  P 

For  fluid  flow  in  microchannel  at  constant  mass  flow  rate,  pressure  gradient  —is  required  to  drive  the  liquid. 
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The  co-efficient  of  flow  resistance  /  is  known  as  friction  factor  and  it  is  related  to  the  Reynolds  number  by  Blasius 
equation  for  turbulent  flow, 

/  =  0.3164(Re)0-25  (5) 

The  Blassius  Equation  gives  the  friction  factor  for  turbulent  flowing  smooth  pipes  for  Re  <  105  [16]. 


For  fully  developed  laminar  flow  at  macroscale  in  pipe  can  be  written  as, 


f  = 


Re 


and, 


(6) 


Re 


pvD 


(7) 


From  conventional  theory,  the  product  of  the  Reynolds  number  and  friction  factor, Re  X  /  =  C  the  value  of  C=64 
for  circular  pipe  flow. 


Table  1:  Friction  Factor  Co-Relations  for  various  Geometrical  Configurations  [17] 


Rectangular 

Friction  factor 

a/b=  1 

56.92/Re 

a/b=2 

62.20/Re 

a/b=3 

68.36/Re 

a/b=4 

72.92/Re 

a/b=6 

78.80/Re 

a/b=8 

82.32/Re 

Triangular 

A=30° 

52.28/Re 

A=60° 

53.32/Re 

A=90° 

52.60/Re 

A=120° 

50.96/Re 

Shah  and  Fondon  correlation  [18]  can  be  applied  for  laminar  flow  for  both  conditions  whether  flow  has  reached 
fully  developed  state  or  not  at  the  outlet. 


A*'--  =  (8) 

Where, 

fRe  =  96(1  -  1.3553a  +  1.9467a2  -  1.7012a3  +  0.9564a4  -  0.2537a5)  (9) 

The  reliability  of  the  above  correlations  has  been  validated  by  the  recent  experimental  work  of  Investigation  of 
liquid  flow  in  microchannels  [19]. 

The  effects  of  friction  factor  and  the  pressure  drop,  for  single  phase  flow  in  microchannels,  over  the  hydraulic 
diameter  range  of  1pm  to  100pm  over  the  mass  flow  rate  range  of  0.5  gram/sec  to  10  gram/sec  are  discussed.  According  to 
the  Blassius  equation  for  friction  factor,  /  =  0.3164Re-0,25  defines  the  dependence  of  friction  factor,  over  the  Reynolds 
number  for  turbulent  flow. 
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Figure  2:  Plot  of  Friction  Factor  vs.  Reynolds  Number  for  Dh  =  lOjum 


Figure  3:  Plot  of  Friction  Factor  vs.  Reynolds  Number  for  Dh  =  40pm 


The  above  figure  2,  shows  the  variation  in  friction  factor,  with  respect  to  Reynolds  number  for  10pm  diameter 
however,  the  figure  3  shows  the  variation  for  40pm  diameter.  The  behaviour  of  friction  factor  characteristic  for  air,  liquid 
nitrogen  and  liquid  hydrogen  is  compared  and  seen  that,  as  the  value  of  Reynolds  number  increases,  the  value  of  friction 
factor  for  all  three  flowing  fluids  decreases.  It  can  be  seen  further  for  both  cases  that  the  liquid  nitrogen  displays  the 
highest  value  of  friction  factor  and  liquid  hydrogen  displays  the  lowest  value  of  friction  factor  among  all  three  fluids. 
One  important  observation  was  also  made  that  as  the  hydraulic  diameter  was  increased  for  constant  mass  flow  rate,  the 
value  of  friction  factor  was  increased.  It  also  can  be  observed  from  the  figures  that  the  variation  in  friction  factor  occurs  as 
the  conventional  channels  on  the  basis  of  moody  chart  and  comes  in  completely  in  agreement  with  conventional  behaviour 
for  friction  factor. 


Figure  4:  Moody's  Diagram 
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Above  given  moody  chart  in  figure  4  depicts  the  behaviour  of  fluid  flow  according  to  the  conventional  theory  and 
it  shows  that  for  turbulent  flow,  the  variation  of  friction  factor  follows  a  curve  in  decreasing  order  as  the  value  of  Reynolds 
number  keep  increasing  and  same  behaviour  is  followed  by  the  air,  liquid  nitrogen  and  liquid  hydrogen  thus  it  can  be 
conclude  on  the  safe  side  that  fluid  flow  behaviour  in  microchannels,  does  follows  the  conventional  theory. 


Pressure  Drop 

Basically,  there  are  three  major  components,  which  have  their  contribution  to  the  overall  pressure  drop  in  fluid 
flow,  in  microchannels  or  duct.  The  losses  at  inlet  and  exit  need  to  be  quantified  as  the  flow  in  microchannels  takes  place. 
While,  evaluating  the  hydrodynamic  developing  length  care  must  be  taken  into  account.  Finally,  the  rest  of  the  length  will 
have  the  fully  developed  frictional  loss.  The  components  of  the  total  pressure  drop,  can  be  seen  in  the  following  given 
equation  [21]. 

Ap  =  [/d  +  ko  +  fapJL]  (10) 


The  co-efficient  of  inlet  and  outlet  loss  are  represented  by  iq  and  k0.  However,  The  Navier-Stokes  (NS)  equation 
can  be  simplified,  in  terms  of  Hagen-Poiseuille  equation  and  the  feature  of  Pressure  drop,  for  the  circular  pipe  flow, 
according  to  the  conventional  theory  can  be  described  as, 


Ap  =  128vLQ 
nDA 


(ID 


Where  D  and  L,  are  the  diameter  and  the  length  of  the  microchannel,  respectively,  AP  is  the  pressure  drop  amidst 
the  two  ends  of  the  microchannel,  and  p  is  the  viscosity  of  the  fluid.  This  equation  is  valid  for  laminar  regime,  as  well  as 
turbulent  regime  for  incompressible  flow. 


Figure  5  :  Plot  against  Pressure  Drop  and  Reynolds  Number,  Dh=  30pm 
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Figure  6:  Plot  against  Pressure  Drop  and  Reynolds  Number,  Dh=  50pm 


In  the  above  figure  5,  the  behaviour  of  pressure  drop,  for  flowing  liquid  air,  liquid  nitrogen  and  liquid  hydrogen 
for  hydraulic  diameter  30  pm  and  in  the  figure  6,  the  behaviour  for  hydraulic  diameter  70pm  is  shown.  From  the  above 
figures,  it  can  be  seen  that,  as  the  value  of  Reynolds  number  for  particular  hydraulic  diameter  is  increased,  with  the 
increasing  value  of  mass  flow  rate,  the  value  of  pressure  drop,  also  keep  on  increasing  for  all  flowing  fluids.  As  a  result,  it 
can  be  predicted  that,  by  pressure  drop  being  the  function  of  mass  flow  rate,  according  to  Bernoulli  equation, 

(Pi-P2)=^[(Z2-Zi)+(X^il+HL  (13) 


It  was  also  observed  that,  as  the  hydraulic  diameter  was  increased  at  constant  mass  flow  rate,  the  value  of  pressure 
drop  decreased  for  all  three  flowing  fluids.  Both  figures  demonstrate  that,  for  every  diameter  air  yields  the  highest  value  of 
pressure  drop  and  liquid  nitrogen  yields,  the  lowest  value. 


In  the  figure  7,  a  graph  shows  a  variation  in  pressure  drop  against  mass  flow  rate  for  different  hydraulic  diameters. 
This  graph  also  shows  that  as  the  value  of  mass  flow  rate  increases,  the  pressure  drop  increases  for  all  considered 
diameters.  However,  it  also  can  be  seen  that  as  the  hydraulic  diameter  increases,  the  pressure  drop  decreases  and  the 
microchannel,  having  hydraulic  diameter  of  1pm  shows  the  least  value  of  pressure  drop,  among  the  considered  diameters. 


Wall  Shear  Stress 


When  the  fluid  moves  along  the  solid  boundary,  it  offers  some  resistance  force,  which  reduces  the  velocity  of  the 
fluid;  however,  the  viscosity  plays  the  role  of  an  important  factor. 

n  =  -k  (13> 

dy 


Figure7:  Plot  against  Pressure  Drop  and  Mass  Flow  Rate 
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No  slip  phenomenon  is  defined,  when  wall  boundary  does  not  offer  any  resistance  to  the  fluid  motion,  thus  fluid 
will  flow  without  friction  at  the  wall  contact. 


Figure  8:  Plot  of  Wall  Shear  Stress  and  Reynolds  Number,  Dh=  30pm 


Figure  9:  Plot  of  Wall  Shear  Stress  and  Reynolds  Number,  Dh=  90 jam 

In  the  above  given  figure  8,  it  shows  the  effect  of  increasing  value  of  Reynolds  number,  over  the  wall  shear  stress 
for  hydraulic  diameter  of  30pm  and  figure  9,  shows  the  graph  for  90pm,  for  the  varying  mass  flow  rate.  Both  figures  show 
that,  as  the  value  of  mass  flow  rate  increases,  Reynolds  number  increases  and  with  the  increasing  value  of  Reynolds 
number,  the  value  of  wall  shear  stress  increases.  It  is  also  observed  that,  as  the  hydraulic  diameter  increases,  wall  shear 
stress  is  decreases.  Both  figures  show  that,  the  air  yields  the  highest  value  of  the  wall  shear  stress  and  liquid  nitrogen  yields 
the  lowest  value  of  wall  shear  stress. 

Heat  Transfer  Characteristics 
Rate  of  Heat  Transfer 

Heat  transfer  by  means  of  convection  process,  provides  a  mode  of  transferring  the  heat  between  fluid  to  fluid  or 
fluid  to  solid,  quickly  from  the  exchange  surface.  Fluids  used  in  the  systems  and  processes;  incur  the  changes  in  their 
thermal  state  in  heat  exchanger.  The  basic  equation  for  heat  transfer  by  convection  is, 

q  =  hA(Ts  -  7»  (14) 

However,  when  a  cold  fluid  is  flowing  in  a  hot  channel  with  a  solid-liquid  temperature  difference,  the  heat 
transfer  between  the  liquid  and  solid  at  boundary  takes  place.  From  conservation  of  energy,  the  rate  of  heat  transfer  imbibe 
by  the  liquid  results  in  a  liquid’s  temperature  increase, ATt 
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q  =  mCpATi 


(15) 


Where,  m  is  the  mass  flow  rate  of  fluid  and  cp  isspecific  heat  of  the  liquid. 

According  to  the  conventional  theory,  the  heat  transfer  co-efficient,  is  unswayed  by  the  Reynolds  number  in  the 
fully  developed  region  of  laminar  flow  is  given  by, 

Nuxk 


h 


(16) 


The  nusselt  number  for  the  flow  in  microchannel,  can  be  calculated  by  Chilton-Colburn  analogy  [22]  given  as, 

Nu  =  0.125 fRe  x  Pr033  (17) 

The  local  convective  heat  transfer  coefficient  and  local  Nusselt  number  [22]  are  expounded  as, 


hx= - 3d - 

N  A ch(Tw,x  -  Tx) 


(18) 


According  to  the  conventional  theory,  the  heat  transfer  characteristics  in  laminar  flow  for  circular  duct,  can  be 
expressed  in  terms  of  nusselt  number,  which  is  Nu  =  Nu  =  4.36  [Ts=  cons.]  and  Nu=  3.66  [q”=  cons.]  [22]. 

However,  for  gas  flow  in  circular  channels,  the  value  of  nusselt  number  depends  upon  the  Reynolds  number,  as 

follows, 


Table  2:  Nusselt  Number  Correlation  for  different  Ranges  of  Reynolds  Number  [22] 


Nusselt  Number 

Reynolds  Number  Range  [Gas  Flow] 

Nu=  0.989Re  °'330Pr  ^ 

Re=0.4-4 

Nu=  0.9 1 1  Re  03s5Pr 0  33 

Re=4-40 

Nu=  0.683Re  °'466Pr 103 

Re=40-4000 

Nu=  0. 1 93Re  06l8Pr 11,3 

Re=4000-40,000 

Table  3:  Nusselt  Number  Co-Relations  for  various  Geometrical  Configurations  for 
Constant  Surface  Temp  and  Constant  Surface  Heat  Flux  [22] 


Rectangular 

Ts=  Constant 

qs=  Constant 

a/b=  1 

Nu=2.96 

Nu=3.61 

a/b=2 

Nu=3.39 

Nu=4.12 

a/b=3 
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Figure  10:  Plot  against  Heat  Transfer  and  Reynolds  Number,  Dh=  20pm 


Figure  11:  Plot  against  Heat  Transfer  and  Reynolds  Number,  Dh=  50pm 

In  the  below  given  figure  10  and  figure  11,  the  graph  for  the  hydraulic  diameters  20pm  and  50pm,  respectively, 
shows  the  effects  of  Reynolds  number  variation,  over  the  surface  heat  transfer.  It  can  be  observed  from  the  given  graphs, 
for  hydraulic  diameters  that,  as  the  value  of  Reynolds  number  increases  with  the  increasing  value  of  mass  flow  rate,  the 
values  of  surface  heat  transfer  increases.  It  was  observed  that,  the  liquid  hydrogen  yield  the  highest  value  of  heat  transfer 
and  air  yield  the  lowest  value  of  heat  transfer,  among  the  all  three  flowing  fluids.  It  also  can  be  seen  for  the  given  both 
graphs  that,  as  the  value  of  hydraulic  diameters  increases,  the  value  of  heat  transfer  also  increases  but  the  increase  in  heat 
transfer  is  not  so  significant.  These  graphs  also  indicate  that,  liquid  hydrogen  and  liquid  nitrogen,  being  the  cryogenic 
fluids  show  astonishing  heat  dissipation  properties,  when  they  are  combined  with  microchannels.  It  also  can  be  hoped  that, 
they  will  be  extremely  useful  tool  in  the  area  of  microelectronics,  for  the  purpose  of  eliminating  the  thermal  restrictions,  in 
order  to  enhance  the  performance  of  electronic  instrument,  as  well  as  in  the  implementation  in  the  area  of  cooling  of 
turbine  blades. 

Nusselt  Number 

For  the  analysis  of  heat  transfer  characteristics,  Nusselt  number  is  the  best  and  most  important  factor,  which  can 
be  taken  into  consideration.  Nusselt  number  is  dimensionless  number,  which  defines  the  dependence  of  heat  transfer  co¬ 
efficient  of  to  the  diameter  of  the  duct  and  heat  transfer  co-efficient,  shows  the  effects  over  the  heat  flux  or  the  rate  of  heat 
transfer. 


Nu  = 


hD 


(19) 
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Figure  12:  Plot  of  Nusselt  Number  and  Reynolds  Number,  Dh=  ljum 
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Figure  13:  Plot  of  Nusselt  Number  and  Reynolds  Number,  Dh=  20  pm 


The  above  given  Figure  1 1  for  the  hydraulic  diameter  lpm  and  figure  12  for  hydraulic  diameter  of  20pm,  explains 
the  effects  of  Reynolds  number,  over  the  Nusselt  number.  It  can  be  observed  from  the  graphs  that,  as  the  Reynolds  number 
value  increases  with  the  increasing  value  of  mass  flow  rate,  the  value  of  Nusselt  number  also  increases.  Figures  show  that, 
by  increasing  the  Reynolds,  the  Nusselt  number  increases  however,  on  comparing  the  both  figures  it  also  can  be  observed 
that,  for  a  particular  value  of  mass  flow  rate  if  hydraulic  diameter  of  the  microchannel  is  increased,  and  the  value  of 
Nusselt  number  decreases.  Both  figures  depicts  that,  the  air  yield  the  highest  value  of  nusselt  number  in  microchannels  and 
the  liquid  nitrogen,  yields  the  lowest  value  for  the  nusselt  number,  for  all  the  considered  diameters. 


Figure  14:  Plot  against  Nusselt  Number  and  Hydraulic  Diameters 
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Figure  15:  Plot  against  Nusselt  Number  and  Hydraulic  Diameters 

Figure  14  is  a  graph  drawn  between  the  nusselt  number  and  the  hydraulic  diameter,  range  from  1pm  to  100pm. 
This  graph  depicts  that,  as  the  diameter  of  the  microchannels  increases,  the  value  of  nusselt  number  continue  to  decrease, 
while  mass  flow  rate  was  kept  constant.  It  also  shows  that,  the  air  yields  the  highest  value  of  nusselt  number  and  the  liquid 
nitrogen,  yields  the  lowest  value  of  nusselt  number,  among  all  three  flowing  fluids.  Same  kind  of  behaviour  can  be  seen 
from  the  figure  15.  However,  on  comparing  the  both  graphs  it  can  be  seen  that,  as  the  mass  flow  rate  is  increased,  the  value 
of  nusselt  number  also  gets  higher. 

CONCLUSIONS 

In  this  paper,  a  computational  investigation  to  understand  the  Thermohydraulic  characteristics  in  microchannels 
by  considering  air,  liquid  nitrogen  and  liquid  hydrogen,  as  flowing  fluid  was  carried  out.  For  single-phase  flow  and 
turbulent  flow  regime,  the  fluid  flow  characteristics,  in  microchannels  with  hydraulic  diameter  range  of  1pm  to  100pm, 
obeys  the  classical  flow  theory.  The  friction  factor  behaviour,  for  all  three  flowing  fluids  at  varying  mass  flow  rate  in  the 
range  of  0.5  gram/sec  to  10  gram/sec,  follows  the  same  as  for  the  conventional  channel,  as  seen  according  to  the  moody 
diagram.  Pressure  drop  was  found  to  be  increasing,  with  increasing  Reynolds  number.  It  was  observed  that,  wall  shear 
stress  increases  with  increasing  value  of  Reynolds  number  and  follows  shear  thinning  behaviour.  In  the  case  of  heat 
transfer  characteristics,  as  the  Reynolds  number  increased  with  increasing  value  of  mass  flow  rate,  the  rate  of  heat  transfer 
was  also  increasing.  It  was  seen  that  nusselt  number,  increasing  with  increasing  Reynolds  number  for  all  three  flowing 
fluids.  It  was  observed  that,  with  the  use  of  cryogenic  fluids  such  as  liquid  nitrogen  and  liquid  hydrogen,  the  heat  removal 
capacity  of  microchannels  can  be  increased  sufficiently,  so  it  could  be  employed  in  the  turbine  blades  and  vanes  as  a 
substitute  of  cooling  technique. 
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